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ABSTRACT 


Submitted  to  the  Department  of  Naval  Architecture  and  Marine 
Engineering  on  May  23,  19&9,  in  partial  fulfillment  of  the 
requirements  for  the  Master  of  Science  Decree  in  Mechanical 
Engineering  and  the  Professional  Decree,  Naval  Engineer. 

Submarine  Pontoon  Supported  Vessels  are  being  used  or  consid- 
ered for  drilling  platforms,  oceanographic  research  ships,  and 
mobile  satellite  tracking  stations  due  to  their  inherently- 
small  notions.   This  thesis  explores  the  possibility  of  reduc- 
ing these  motions  still  further  by  using  active  displacement 
control.   The  equations  of  motion  for  a  typical  model  are 
derived  and  solved  for  heave,  pitch,  and  roll.   The  forces 
required  to  reduce  the  motions  are  optimized  and  the  result- 
ing motions  are  determined.   .he  power  required  to  maintain 
a  specified  motion  amplitude  is  derived. 

The  reduction  of  heave  amplitude  of  a  vessel  near  resonance 
requires  excessive  power,  indicating  that  the  vessel  should 
be  designed  with  a  natural  period  above  that  of  the  expected 
waves.   However,  the  use  of  active  displacement  control  to 
provide  essentially  zero  heave  amplitude  in  a  moderate  sea 
state  is  potentially  attractive  due  to  the  small  average  power 
requirement. 


Thesis  Supervisor:     Forbes  T.  Brown 

Title:   Associate  Professor  of  Mechanical  Engineering 


-  2  - 


TABLE  OF  CONLLLI'S 

PAGE 

TITLE  PAGE  1 

ABSTRACT  2 

TA3LE  OF  CONTENTS  3 

NOTATION  k 

I  INTRODUCTION  6 

II  DESCRIPTION  OF  MODEL  10 

III  DERIVATION  OF  FORCES  AND  EQUATIONS  OF  iiCTION  13 

IV  SOLUTION  OF  EQUATIONS  OF  MOTION  21 

V  OPTIMAL  CONTROL  22 

VI  POWER  REQUIREMENTS  25 

VII  DISCUSSION  AND  CONCLUSIONS  29 
APPENDIX  A  DETAILS  OF  FORCE  DERIVATION  3^ 
APPENDIX  B  DETAILS  OF  EQUATIONS  OF  OTICN  37 
AxPENDIX  C  DETAILS  OF  SOLUTION  *fl 
APPENDIX  D  DETAILS  OF  CONTROLLED  MOTION  hk 
APPENDIX  E  DETAILS  OF  POWER  CALCULATIONS  V? 
REFERENCES  51 


-  3  - 


:',<    r.  i'.OII 

A  Incident  wave  amplitude 

b  Eleight  of  the  SPSV's  center  of  gravity  above  calm 

waterline 

c  subscript  applying  to  a  column 

C  Drag  coefficient  of  a  column 

c  Added  mass  coefficient 
am 

c ,  Drag  coefficient 

d  Diameter  of  a  column 
c 

d,  Diameter  of  a  hall 
h 

H  Wetted  length  of  a  column 

h  Subscript  applying  to  a  hull 

i  Subscript  referring;  to  a  particular  segment 

I  ,1  ,1  Platform  moments  of  inertia  about  x,  y,  z   axes 

k  V/ave  number  =  2  /L   =  w  /g 

w 

L  Length  of  a  hull 

L.  Length  of  a  hull  segment  =  \/2 

L  vave  length 

M  3PSV  mass 

m  displaced  column  mass 
c  ^ 

m.  displaced  hull  mass 

Am  Added  mass  or  hvdrodynamic  inertia  =c 

am 

S  cos(kx      cos*  )    cos    (ley      siruX  ) 
c  c 

t  time  variable 

T  period  =  1/w 

X,YfZ  Force  components 

x,y,z  Cartesian  coordinate  system 

x  . -  y  .   .ocation  of  axes  of  i-th  coordinate  system 
ci'  ■'ci 

z,  vertical  distance  of  a  hull  axis  below  the  SPoV's 
n 

center  of  gravity 

<*  Incident  wave  angle 

£  Damping  parameter 
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V  «  > y  Linear  displacement  variables 
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I.   INTRODUCTION 

In  recent  years  a  considerable  amount  of  attention 
has  been  focused  on  floating  platforms  for  a  wide  variety 
of  missions.   These  platforms  have  been  identified  by 
many  -names,  including  "Mobile  Column  Stabilized  Platforms 
(MCSP)"  and  "Stable  Ocean  Platforms  (SOP)".   For  purposes 
of  this  paper  they  will  be  referred  to  as  "Submarine 
Pontoon  Supported  Vessels  (SPSV)".   The  SP3V  is  charac- 
terized by  its  relatively  small  wave  induced  motions. 
It  consists  of  a  platform  type  structure  supported  by 
columns  attached  to  submarine  hulls.   The  main  buoyancy 
elements  are  the  hulls,  and  since  they  are  located  below 
the  area  of  major  wave  effects,  the  motions  of  the  vessel 
are  small.   The  supporting  columns,  which  are  exposed  to 
wave  action,  are  widely  spaced,  small  in  relation  to  the 
wave  lengths,  and  present  the  same  minimum  cross  section 
and  waterplane  area  to  the  sea  in  any  orientation.   As  a 
result  of  these  characteristics,  the  SPSV  is  useful  for 
missions  which  require  a  ship  to  operate  at  a  specific 
point  in  the  ocean  with  minimum  motions.   Speed  in  trans- 
it may  be  important  but  remains  secondary  to  the  motion- 
free  characteristics.   Since  the  major  working  areas  are 
not  in  direct  contact  with  the  surface  of  the  sea,  the 
design  and  utilization  of  the  upper  structure  need  not 
be  restricted  by  hydrodynamic  considerat i one,  as  is  the 
case  with  a  conventional  surface  ship.   Thus  unconven- 
tional arrangements  and  unique  configurations  are  pos- 
sible to  fulfill  specific  tasks. 

The  oil  industry  has  been  the  first  to  design, 
construct,  and  operate  a  SPSV.   Deep  water  drillin 
requires  a  floating  structure,  and  the  SPSV  has  filled 
this  need.   The  platform  j  rovides  sufficient  space  for 
the  drilling  rig  and  associated  equipment,  and  the  small 
motion  of  the  vessel  permits  drillin   in  low  sea  states 


without  problem.   Most  of  the  existing  platforms  are  non-self 
propelled,  and  they  are  positioned  at  sea  by  anchoring  systems. 
However,  as  the  need  for  drilling  in  deeper  water  increases, 
the  platforms  will  be  equipped  with  dynamic  positioning  devices 
to  maintain  station. 

The  use  of  the  SFSV  for  scientific  research  holds  great 
potential.   The  platform  provides  considerable  volume  and  deck 
area  for  specialized  equipment  and  scientific  facilities.   Due 
to  the  considerably  reduced  motions  as  compare:;  with  a  conven- 
tional ship,  personnel  comfort,  and  hence  effectiveness,  is 
high.   As  a  further  consequence  of  reduced  motions,  laboratory 
facilities  can  be  more  sophisticated,  and  personnel  can  contin- 
ue to  perform  research  in  weather  conditions  which  would  pro- 
hibit work  in  a  conventional  ship.   The  i'C'ICLE  vessel,  although 
not  actually  constructed,  underwent  considerable  design  effort. 
The  vessel  was  intended  to  drill  through  the  earth's  mantle  in 
very  deep  water.   The  vessel  consisted  of  a  large  rectangular 
platform  supported  on  two  submerged  hulls  by  six  cylindrical 
columns.   The  design  called  for  a  propulsion  system  capable 
of  approximately  10  knots,  and  a  positioning  system  using 
ducted  thrusters  in  the  hulls. 

The  National  Aeronautics  and  Space  Administration  (NASA) 
has  shown  interest  in  a  SPSV  for  use  in  missile  and  space 
vehicle  tracking  and  communications.   The  requirements  for 
ballistic  missile  and  space  missions  are  strict,  and  since 
these  missions  often  require  telemetry  and  tracking  coverage 
over  ocean  areas,  a  sea-going  tracking  facility  is  necessary. 
Present  range  tracking  ships  are  significantly  limited  in 
that  the  pitching  and  rolling  of  the  ship  create  radar  data 
inaccuracies  for  which  adequate  compensation  cannot  be  provided. 
In  order  to  minimize  the  effects  of  wind,  current,  and  waves, 
present  tracking  ships  must  normally  maintain  headway,  cover- 
ing a  path  which  further  complicates  the  tracking  problem. 
Crew  discomfort  during  adverse  weather,  and  lack  of  space  for 
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recreation  become  significant  problems  during  extended  missions. 

The  SPSV  overcomes  these  disadvantages  and  provides  other 
unique  advantages.   Since  the  motions  of  the  SPSV  are  consid- 
erably less  than  those  of  present  ships,  data  accuracy  is 
improved,  equipment  reliability  is  increased,  and  personnel 
discomfort  is  alleviated.   The  above  water  platform  can  be 
optimized  to  suit  the  mission,  and  radars  and  other  electronic 
equipment  can  be  located  more  effectively  than  on  a  conven- 
tional ship  form.   Additional  space  is  also  available  for 
living  areas  and  recreation.   The  SFSV  would  not  be  required 
to  steam  in  a  pattern  because,  since  it  is  relatively  insensi- 
tive to  v/ind  and  wave  direction,  a  dynamic  positioning  system 
would  permit  "hovering"  over  a  specific  spot.   ihis  capability 
implies  that  the  endurance  of  the  SPSV  on  station  would  be 
greater  than  that  of  a  conventional  ran^e  tracking  ship. 

The  SPSV  has  potential  as  a  support  vessel  for  deep 
diving  submersibles,  either  for  oceanographic  research  or  for 
submarine  rescue.   Small,  deep  diving  submersibles  have 
limited  range  and  require  close  support.   The  3, SV  provides 
a  relatively  motionless  vessel  which  could  house  or  support 
such  submersibles. 

There  are  a  number  of  other  possible  missions  for  the 
SPSV,  although  some  may  be  well  in  the  future.   The  salvage 
potential  is  large  due  to  the  capability  of  lifting,  large 
weights  without  developing  large  pitch  or  roll  moments. 
Underwater  structures,  habitats,  and  arrays  will  require 
substantial  tending,  particularly  during  construction,  and 
the  SPSV  has  a  unique  capability  to  provide  such  support. 
Further  uses  such  as  a  helicopter  landing  platform  or  a 
sea  going  hotel  are  conceivable. 

With  all  of  the  fore  voing  missions,  the  ability  of  the 
SPSV  to  be  located  at  a  fixed  point,  with  small  motion,  and 
yet  be  capable  of  moving  from  one  location  to  another  are 
primary  considerations.   In  many  cases  it  might  be  desirable 


to  minimize  notions  ctill  further.   This  aight  be  necessary 
to  provide  an  essentially  r.-.otionless  platform,  or  to  extend 
the  working  ran'-e  of  the  SI  3V  into  higher  state  seas.   As 
the  natural  period  of  heave,  pitch,  and  roll  are  large  for 
the  S:  3\T ,    it  might  be  feasible  to  control  Motions  with  the 
use  of  an  active  displacement  control  system.   i'he  system 
envisioned  would  consist  of  additional  ballast  tanks  and 
the  pumping  capacity  to  offset  wave  induced  motions  by 
alternately  increasing  and  decreasing  ballast.   If  the  water 
were  added  and  removed  at  the  proper  location  and  with  the 
proper  phase  relation  to  the  exciting  wave  forces,  the 
motions  of  the  SL SV   in  heave,  roll,  and  pitch  could  theo- 
retically be  reduced  to  near  zero. 

This  investigation  considers  the  development  of  a 
active  displacement  control  system  for  a  SF3V.   To  determine 
the  feasibility  of  such  a  system,  the  power  requirements  are 
determined  for  three  situations.   The  first  is  the  mainte- 
nance of  essentially  zero  motions  in  low  to  moderate  seas. 
The  second  is  the  significant  reduction  of  heave  in  neavy  seas, 
Finally  the  reduction  of  heave,  at  resonance,  to  a  reasonable 
level  is  considered.   The  alteration  of  the  il'SV's  natural 
period  is  also  discussed. 


II.  DESCRIPTION  OF  SPSV  MODEL 

The  SPSV  consists  of  two  cylindrical  hulls  which  are 
joined  to  the  working  platform  by  four  cylindrical  columns. 
The  working  platform  is  rectangular  in  shape.   Figure  1  shows 
a  simplified  lines  drawing  of  the  SPSV. 

A  Cartesian  coordinate  system  (x,  y,  z)  is  established 
with  the  origin  located  at  the  undisturbed  center  of  gravity 
(KG).   The  x-axis  is  positive  towards  the  bow,  the  y-axis  is 
positive  to  port,  and  the  z-axis  is  positive  upwards.   The 
undisturbed  water  surface  is  the  plane  z=  -b. 

It  is  assumed  (ref.  1)  that  the  columns  are  far  enough 
apart  so  that  each  column  and  the  portion  of  the  hull  directly 
beneath  can  be  considered  a  hydrodynamically  separate  segment, 
and  any  interactions  between  segments  can  be  neglected.   The 
segments  are  numbered  as  shown  in  figure  2. 

The  hull  axes  are  parallel  to  the  x-axis  and  the  column 
axes  are  parallel  to  the  z-axis.   For  convenience,  an  addi- 
tional coordinate  system  (x. ,  y. ,  z.)  is  introduced  in  the 

i    segment  with  the  origin  at  (x  . ,  y  . ,  0)  of  the  primary 

ci    ci* 

coordinates. 

Displacements  of  the  center  of  gravity  are  described  by 
surge  (  \   ),  sway  (  tj  ) ,  and  heave  (  \  ),  and  rotations  about  the 
x,  y,  and  z  axes  by  roll  ($),  pitch  (  <9 ) ,  and  yaw  ( v[>  ) . 
Displacements  of  a  point  located  on  a  column  axis  (x  . ,  y  . ,  z. ) 

r  CI    CI    1 

can  be  related  to  displacements  of  the  center  of  gravity  by: 

li"  \    +  Zi9"  :'ci^  (1"a) 

\i   -    V  xoi^  -  \$  (1-b; 

Si-  S  +  y0i^-xci6  (1-c) 
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r^\ 


r^ 


r~\ 


v_y 


Figure   1      SPSV   Lines 
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I  * 


ir;ure   2     Segment    Numbering 
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III.  DERIVATION  OF  FORCES  AND  EQUATIONS  OF  MOTION 

Wang  (ref.  1)  has  derived  the  motions  and  forces  of  a 
floating  platform.   The  derivation  of  the  motions  and  forces 
of  the  SPSV  draw  heavily  from  this  work. 

It  is  assumed  that  the  SFSV  is  subjected  to  a  regular 
train  of  two  dimensional  deep  water  gravity  waves  whose 
velocity  potential  is  of  the  form: 

T    ,.  w  k  ( z  +  b  )  .   ~  .  _  v 

<p  -   -A~  e        sin  X  (2) 

From  che    theory  of  gravity  waves  (ref.  3i  *0 1  the  wave 

profile  can  be  written: 

^   1  c^  j»    .  k(z  +  b)    r  r_v 

<o  =  -      =  Ae        cos  L  (3) 

8    b  t 

where : 

L  =  k(x  cos  Cn.  -  y  sin  cO  -wt 

g  =  gravitational  constant 

t  =  time  variable 

A  =  surface  wave  amplitude 

2 
k   =  v/ave   number    =   2  Tf/L      =  w  /g 

w      ° 

L  =  wave  length 
w 

w  =  circular  frequency  of  the  wave 

ck   =  angle  of  the  wave  front  measured  clockwise  from 

o 
the  negative  x-axis.   (For  port  beam  seas,  <*  =  90  , 

for  head  seas,  c<  =  loO  ). 

It  is  furlher  assumed  that  the  column  and  hull  diameters 
are  small  compared  to  the  wave  length,  therefore  the  motion  of 
the  water  across  the  diameters  can  be  considered  uniform. 

The  forces  acting  on  the  SPSV  are  hydrostatic  and 
hydrodynamic.   The  hydrostatic  force  is  the  restoring  force 
expressed  as: 

Hydrostatic  Force  =  -r   fn  d  (AcoslT'.  -  Si ) 
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The  hydro dynamic  forces  are  c  msed  by  changes  in 
pressure  on  the  segment  due  to  relative  velocities  and 
accelerations.   ..hose  hydrodynamic  forces  in  phase  with 
accelerations  are  termed  intertial,  and  those  in  v.hase  with 
velocities  are  termed  damping. 

The  inertial  forces  consist  of  a  wave  force  and  a  force 
due  to  the  added  mass  of  the  segment.   The  wave  force  is  due 
to  the  pressure  on  the  segment  resulting  from  the  segment's 
interference  with  the  acceleration  of  the  water  particles. 
Korvin  -  Kroukovsky  (ref.  *+,  cnapter  2)  expresses  this  force 
as: 

Wave  Force  =  -(1  +  k)  (pressure  gradient)  (body  volume)   (5) 


From  gravity  wave  theory,  the  pressure  is: 


?    = 


1ft   =    *9 


(6) 


Substituting  equation  (3)  into  equation  (6),  the  expression 
for  pressure  becomes: 

ft.k(z+b)/—  (•?) 

p  =  ^>g*e       "cos  X.  Kn 

The  pressure  gradient  is  easily  calculated  from  equation 

(7). 

In  equation  (5),  k  is  the  coefficient  of  accession  to 
inertia,  and  1  +  k  is  defined  in  reference  k%    chapter  2  as: 

total  inertia  of  a  bod;/  floating  in  a  fluid     ,  g  x 
inertia  of  the  fluid  displaced  by  the  body 

Equation  (8)  can  also  be  written: 

1  +  k  =  ^— + — -   ,  where  ftm  is  the  added  mass  or  hydro- 
m 

dynamic  inertia. 
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The  inertia  force  due  to  the  hydrodynaraic  inertia  is  the 
product  of  the  hydrodynaxnic  inertia  and  the  relative  accelera- 
tion.  The  damping  force  is  the  product  of  the  damping 
coefficient  and  the  relative  velocity  between  the  water  and 
the  segment. 

The  water  velocities  in  the  x,  y,  and  z  directions  can 
be  obtained  from  the  velocity  potential,  equation  (2),  as 
follows: 

V  =  -  — —  =  Awek(z  +  b)cos<*cosT  (9-a) 

c^  X 

V  =  -  A£   =  -Awek(z  +  b)sin<*  cost  (9-b) 

y        by 

c>  j>    .   k(z  +  b)    _.  f-      x 

V  =  -  —  =  Awe        sin  X  (^-c) 

z      2>L 

The  water  accelerations  can  be  determined  by  taking  the  time 
derivatives  of  the  velocities,  and  are: 

,  2  k  ( z  +  b )    _  .  _  / -.  r       \ 

V  =  Aw  e         cos«.sxnt  (10-a) 
x 

V  =  -Aw2ek(z  +  b)sino<  si-iT  (10-b) 

y 

w      .  2  k  ( z  +  b  )    ~.  / , n  \ 

V  =  -Aw  e        cos  <-  (,10-c; 
z 

The  relative  velocities  and  accelerations  can  be 
obtained  by  calculating  the  water  velocities  and  accelerations 
at  1;r.e  desired  segment,  and  then  subtracting  the  segment 
motions . 

The  forces  acting  on  each  segment  can  be  resolved  into 
six  components  : 

a   is  the  force  on  the  i-th  column  in  the  x-direction. 
ci 

V  .  is  the  force  on  the  i-th  column  in  the  y-direction. 
ci 

Z   is  the  force  on  the  i-th  column  in  the  z-direction. 

ci 
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X,  .  is  the  force  on  the  i-th  hull  in  the  x-direction. 

Y,  .  is  the  force  on  the  i-th  hull  in  i:he  y-direction. 
hi  J 

Z,  .  is  the  force  on  the  i-th  hull  in  the  z-direction. 
hi 

The  hydro dynamic  inertia  for  the  motion  of  a  cylinder 
along  its  longitudinal  axis  is  negligible.   The  hydrodynamic 
inertia  for  the  transverse  motion  of  a  slender  vertical 
cylinder  can  be  approximated  by  its  displaced  mass  (ref.  1). 
The  hydrodynamic  inertia  for  the  transverse  motion  of  a 
horizontal  cylinder  below  a  free  surface  is  a  function  of 
the  oscillation  frequency  and  submergence  ratio  (depth/radius) 
of  the  cylinder,  and  is  plotted  in  figure  6  of  reference  1. 

The  damping  force  is  composed  of  wave  damping  and  eddy 
damping.   The  wave  damping  is  associated  with  progressive 
waves  radiating  outward  from  the  body,  and  it  is  linear  with 
respect  to  the  relative  velocity.   The  wave  damping  force  due 
to  the  columns  is  negligible  and  a  wave  damping  parameter  for 
the  hull  is  plotted  in  figure  7  of  reference  1  as  a  function  of 
oscillation  frequency  and  submergence  ratio.   The  wave  damping 
coefficent  is  then  expressed  as: 

C   =  S   m,  w  (11) 

w      h 

The  eddy  drag  is  associated  with  the  generation  of 
turbulent  eddies  around  the  segments  and  is  nearly  proportional 
to  relative  velocity  squared.   The  motion  responses  of  the 
platform  are  small  in  the  range  outside  of  resonance,  therefore 
this  damping  is  small  and  can  be  neglected,  however  near 
resonance  it  must  be  included.   The  eddy  drag  normal  to  a 
cylinder  in  cross  flow  is: 

D  =  cj>  f  o\V|V  (12) 

where:   c,  =  drag  coefficient 
d 
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o     =   water  density 

j   =  projected  area  of  the  cylinder 

V  =  relative  velocity 

The  absolute  value  of  V  is  used  in  equation  (12)  to  show 
that  the  force  is  always  in  the  direction  of  the  relative 
velocity.   The  value  of  c,  is  determined  experimentally , 
and  an  average  value  of  1.0  seems  reasonable  (ref.  1). 

The  relative  velocity  in  equation  (12)  can  be  written 
in  the  form: 

V  =  w  Xcos(wt-  £  )  (13) 
where : 

w  =  frequency  of  oscillation 
X.  =  amplitude  of  oscillation 

£  =  jjhase  angle 
oince  the  velocity  is  periodic,  the  eddy  drag  can  be  linear- 
ized using  the  describing  function  technique  (ref.  5,  6), 

and  written  in  the  form:  D   =  C  V.   The  describing  function 

v    v 

approach  assumes  a  periodic  input  (V)  and  further  assumes 
that  the  only  significant  component  of  the  output  (drap;) 
is  the  component  at  the  in^  ut  frequency. 

The  describing  function  can  be  written  as  the  ratio 
of  the  fundamental  drag  component  to  the  velocity  amplitude. 

C  =  -i  (14) 

v   wx 

The  fundamental  draft  component  can  be  determined  from  the 
Fourier  series  expansion  of  equation  (12).  The  resulting 
expression  for  the  linearized  eddy  drag-  coefficient  is: 

c      =  St?  cJiPSw  X  (15) 

v        3TT     d    » 

The  linearized  damping  coefficients  for  the  i-th  seg- 
ment can  then  be  written  as: 

-  17  - 


L 
0    .    =  ■=—-  c,Pd    liwX  (16-a) 

L 

yi        37T      d»    c       'cyx  (l6-b) 

H    .    =Sm,w   +  ~^c,-P  d,    J&wX,    •  (16-c) 

yi  h  3TT    d  i     h  '  hyx 

r  4 

zi       °h  3  TT     d  \    h  'lizx  (16-d) 


"Y.       .    is    the   mean   amplitude    of   displacement    of   column   i 

along  the   x-direction.     V        .,   V.,     .,    V-i     .can   be   interpreted 

to  ^  cyx'  r    hyx'    hzi  ^ 

in  a  similar  manner. 

The  details  of  the  calculation  of  force  components, 
as  adapted  from  reference  1,  are  shown  in  appendix  A,  and 
are  expressed  as  follows: 


X  .  =  -(m  k  .  +  J  .  \      .  )  +  3m  w  A  j   cos  <*  sin  T  . 
ex       c  i  ci     xi  »  ci       CO  1 

+  C  .wAQ  coscXcost.  (17-a) 

xi    o  1 

Y   =  -(m  y\      .  +  C  .  "H  .  )  -  3m  w  A 0  sin^sint. 
ci      c  l ci    yi  *  ex      c     o  x 

-C    .wAA    sinotcosT.  (17-b) 

yi       -o  x 

2ci  =  -  ¥?**»  I±  -  'v^o  "  W  >  AcosTi  (17-c) 

X.  .    =   m,  w   A  Ucosc(sinX  .  (17-d) 

hx  h  I  x 

"*  2 

Y,  .    =   -(Am,  Yl  ,  .    +    -    .  T\,  .  )    -    (m,     +   2  Am,  /w   AusindsinT 
hx  h  Ihx  yx    Lhi  h  h  »  x 

-    H    .  wAUsin  ok  cos  T  .  (17-e) 

yx       T  x 

2 
Z.  .    =   -(Am,   <  .     +    H       T     )    "    (mv,  +   2  Am,  )w   AUcosT. 
hx  h0x  zi  Ji  h  h  P  x 

+   H    .wAUsint.  (17-f) 

Zl         "  X 
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."he  equations  of  motion  are  determined  by  summing  the 
forces  and  moments  for  eac  .  segment.   ihe  moments  ^re 
determined  by  calculating  the  force  per  unit  length,  multi- 
plying by  the  moment  arm,  and  then  integrating  over  the 
surface,  being  careful  to  include  buoyancy  effects  in  the 
moment  calculations  (ref.  1). 

4- 
M$    -    £>Cci    +   Xhi>  (l8"a) 

MH    "    t(Yci    +   Yhi}  (l8-b) 

Mi    =1  (2C.    +   Zhj,)  (18-c) 

^  =  Crt"    S(  ^  "ci)z  dz  +  Yhi'h  +  (Zci  +  Zhi)7ci 

+  (nhizh  -  5?^^  dz)c^j  (ib-d) 

I  0=  >  \    \(  4-  x  •  )*  dz  -  X  .z.  -  Z  .x  .  - 

y    £ —  I  £   dz   ex         ci  h    ci  cx 

J  (  7-  Z.  .  )x  dx  +  (m,.z.  -  S?  f  d2.z  dz)0{  (1  -e) 
*L.   dx   nx  hx  h    J  *f  A  ex       J 


I 

z 


4>  =  T  J"  -(X  .  +  X,  .  )y  .  +  Y  .x  .  +  \   (  —  Y   )x  dxi 

T   £-_  J    cx    ai   cx    cx  cx    )  dx   hx      J 


-f) 


..'hen  equation :.  (17)  are  substituted  into  equations  (13),  the 
equations  of  motion  become: 

all  \    +  a15  ©  *  bll  X   +  b15  *  =  J'l  (1  "a) 

"22^  +a2',  f  +  "22  V  b2**t  =  F2  (19"b) 

a33  S  *  b33  5  +  =  33  3  -  h  (1>C) 
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a24  n  +  %k  f +  hzk  \  +  v  <p +  c^  «?  = F. 


a15  i    +   a55 


9+  b15  £  +  b.5  9+  c55  0  =  F, 


+  +  V,^  H>  = 


l66  T  -  "66 


(l.-d) 
(19-e) 
(19-f) 


where  the  coefficients  are  defined  in  appendix  B, 
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IV      SOLUTION   0?    JJiJAriCNS    OF   i  .C  . ;  j    . 

Since  the  control  of  heave,  roll,  and  pitch  are  of 
primary  interest,  the  Motions  of  the  3PSV  in  surge,  sway, 
and  yaw  can  be  eliminated  from  equations  (19)i  resulting 
in  three  uncoupled  differential  equations.   As  outlined 
in  appendix  B,  these  equations  are  of  the  form: 

•»       • 

a,,  *T  +  b-,.,  T.  +  c,,  \  =  f,  sin  wt  +  f_  cos  wt  (2t»-a) 
33  J    33  *    33  *    j>s         3c 

•• •     ■  «•      , 

C   9+  C  2  Q+  C   ©+  C  Q  =   J.  sin  wt  i  J   cos  wt  (20-b) 

C,3$*  °22  f +  C21  f+  C20<P=  J2sSin  Wt  +  J2cCOS  Wt  <2°-c) 

Equations  (20)  can  be  solved  to  determine  the  steady- 
state,  uncontrolled  motions.   Details  of  the  solution  are 
given  in  appendix  G.   For  a  given  frequency,  the  motions 
of  the  ol-SV   are: 


^  =  "S   cos(wt-  £     )  (21-a) 

9=    9      cos(wt-  ia  )  (21-b) 

o  w 

<D  =    (Dq   cos(wt-£^>)  (21-c) 

"^  ,  0  ,  and  <D   are  the  heave,  pitch,  and  roll  amplitudes, 
and  £j,  le  ,  and   £<*>  are  tue  corresponding  phase  angles 
relative  to  the  wave. 
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V.  OPTIMAL  CONTROL 

The  desired  control  may  be  specified  in  many  ways. 
One  is  simply  a  specification  of  the  maximum  allowable 
motions.   This  method  has  the  advantage  of  simplicity  in 
that  an  upper  limit  on  the  motion  is  specified,  and  the 
designer  is  faced  with  a  specific  criterion.   However,  the 
advantage  of  simplicity  may  become  a  disadvantage  since  the 
specification  of  a  certain  motion  may  be  too  expensive  in 
terms  of  providing  the  required  control.   Therefore,  rather 
than  specifying  an  upper  limit  on  motion,  it  is  desirable 
to  determine  an  optimal  control. 

One  method  of  accomplishing  this  is  to  specify  a 
performance  index  in  terms  of  motion  and  the  control, 
minimize  the  performance  index  with  respect  to  the  control, 
determine  the  optimal  control,  and  then  determine  the 
resulting  motions.   One  useful  performance  index  is  based 
on  a  cost  function  which  is  the  sum  of  control  cost  and 
motion  cost. 

The  cost  of  control  can  be  expressed  in  terms  of  the 
capital  cost  required  for  installation  of  the  pumping 
facilities,  etc.,  the  cost  of  operation,  or  some  weighted 
combination  of  the  two.   Since  the  factors  involved  in  the 
determination  of  cost  of  control  can  be  estimated  fairly 
accurately,  an  analytic  expression  for  cost  of  control,  as 
a  function  of  control  force,  can  be  determined. 

The  cost  of  motion  is  less  well  defined  and  depends 
strongly  on  the  mission.   In  determining  this  cost,  the 
owner  must  decide  how  a  particular  motion  will  effect  his 
operation  (drilling,  tracking,  etc.).   Then  he  must  trans- 
late the  effect  into  a  cost  indicator,  expressed  as  a  function 
of  motion. 

The  total  cost  function  can  be  expressed  in  the 
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form: 


C  =  f,  (motion)  +  f?  (control  Action) 

often  the  two  functions  can  be  simplified,  and  the 

cost  function  becomes: 


C  =  qxX  +  pxU  (22) 


Where  q   is  the  cost  of  motion  and  p   is  the  cost  of  control. 
U  is  the  control  force  which  is  to  be  optimized.   X  is  the 
motion  of  the  SPSV,  and  can  be  a  displacement  from  an 
equilibrium  position,  a  velocity,  an  acceleration,  or  any 
combination. 

A  particular  performance  index  is  a  quadratic  form, 
integrated  over  one  period. 
T 

p.i.  =    3[(qx*)2  +  (pxu)2]  dt  (23) 

o 

It  is  further  assumed  that  the  control  force  may  be  expressed 
analytically  in  the  form: 

U  =  U   cos  (wt  -  S  ) 
x  "x 

The  motions  are  also  of  this  form,  X  =  X  cos  (wt  -  £  ) 

c  xc 

where  the  subscript  c  refers  to  the  fact  that  these  motions 
are  those  occurring  with  the  control  forces  applied. 

The  performance  index  can  be  integrated  and  expressed 
in  terms  of  the  notion  amplitude,  control  force  and  costs. 
The  result  of  integrating  equation  (23)  is: 

*•*•  -  t  (<*2  x/  ♦  ^  «x2)  «*> 

To  minimize  the  performance  index,  the  derivative  of 
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equation  (2k)   with  respect  to  t  e  control  is  equated  to 
zero,  giving: 

2TT   2  .   dX    2TT  2n     _  ,._. 

q   a   -777-  +  c   U   =  0  (25) 

w   x  c  dU     w  - x  x 

X 

Equation  (25)  can  be  solved  for  the  optimal  control  force, 

U   ,  with  the  boundary  condition  that  v/hen  the  control 

force  is  zero,  the  amplitude  of  the  motion  is  X   (equations 

o 
(21)).   The  optimal  control  force  is  then: 

U   =  —  (X2  -  X2)^  (26) 

x   p    o    c 
*x 

Since  U  is  an  additional  force  applied  to  the  SPSV 
along  with  the  wave  force,  it  can  be  added  to  the  right 
hand  side  of  the  equations  of  motion.   Equations  (20)  can 
be  solved  to  determine  the  controlled  motions  of  the  SPSV. 
The  details  of  the  calculation  are  presented  in  appendix 
D.   The  result  is: 

A2  -  K  (1  -  A2)^  cos(tf-  &)  -  1  =  0  (27) 


where : 

/\  =  amplitude  factor  =  ratio  of  controlled  to  un- 
controlled motion. 
K   =  a  factor  which  includes  the  ratio  of  cost  of  motion 
to  cost  of  control,  frequency,  and  physical 
characteristics  of  the  3PSV. 
6   =   phase  angle  of  the  control  force 
V=  phase  an^le  of  the  wave  force. 

Equation  (27)  gives  the  motions  of  the  SP3V  resulting  from 
t  e  application  of  an  optimal  control  force  based  on  partic- 
ular motion  and  control  costs.   This  motion,  in  the  form  of  an 
amplitude  factor,  is  used  to  determine  the  power  requirements. 

-  24  - 


vi  ro'/izvi   ?.i:'-„i,']Ha-iEirrs 

The  power  for  the  control  system  must  provide  for  the 
pulping  of  enough  water  at  a  rate  which  will  maintain  the 
motions  of  the  3P3V  within  the  specified  limits.   The 
optimization  of  the  control  force  has  resulted  in  maximum 
allowable  controlled  motions  in  heave,  pitch,  and  roll. 

To  remain  within  the  heave  specification,  added  inertia 
must  be  provided  to  compensate  for  the  uncontrolled  heave 
in  excess  of  that  specified.   This  amount  of  inertia  can 
be  expressed  as: 

Tr  2 

Where  -f  is  the  water  density,  r   d   is  the  cross-sectional 

area  of  one  column,  X   -  \    is  the  excessive  heave 

'   ao    ->c 

which  must  be  controlled,  and  the  expression  is  multiplied 
by  k   to  include  all  four  columns  of  the  SPSV.   The  bar 
over  the  symbol  for  motion  indicates  a  time-varying,  per- 
iodic function,  while  the  absence  of  a  bar  denotes  the 
magnitude,  (e.g.  "S  =  "S   cos(wt  -  i.  ^   )  ). 

To  remain  within  the  roll  and  pitch  specifications, 
added  moment  must  be  provided  in  a  similar  fashion. 
The  amount  of  roll  moment  to  be  added  is: 


I-i 


rr  ,2  2 


<*> 


-*f  i<<  (*o-  *0>  (29) 


The  amount  of  pitch  moment  to  be  added  is: 

Ha  =  k  f   £  d2  x2  (  G      -    9    )  (30) 

9  «  M-   c   c     o     c 

The  required  power  is  the  product  of  the  required  mass 
flow  rate,  and  the  head  against  which  the  pumps  operate. 
From  equations  (28),  (29),  and  (30),  the  mass  flow  rate  in 
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heave,  roll,  and  pitch  are: 


The  head  at  an  arbitrary  point  a  distance  z   feet  below  the 
still  water  line  outside  of  the  hulls  is: 


II      =   z      +   Ae"kz°cos  T-"S       +  x  0     -  y  <f>  (34) 

The  second  term  in  equation  (3^)  includes  the  effect  of  the 
waves,  and  the  last  three  terras  include  the  motion  of  the 
SPSV. 

The  head  inside  of  the  tank  is: 

Pm 

H.  =  ~  +  z.  (35) 

i   fg    i 

z.  is  the  water  height  inside  the  ballast  tank  above  the 
pump,  and  p^  is  the  air  gauge  pressure  in  the  tank.   The 
height  of  the  water  in  the  tank  can  be  expressed  in  terms 
of  the  mass  of  water  in  the  tank  as: 


where  A^  is  the  tank  cross  section  area  and  is  assumed  not 

to  vary  with  z. .   m  is  the  mass  of  the  water  in  the  tank 

at  any  time,  as  determined  from  equations  (28) ,  (29),  or  (30) 

The  total  pump  head  is  the  difference  between  the 
external  and  internal  heads. 

tt  „  -kz„    /v    ~ZT  n      J)  T    m 

H_  =  z.  +  Ae    "cos  t   -  X   +  x  0   -  y  <P -  j-  -  jr 


(37) 
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In  heave,  water  must  be  added  (or  removed)  so  as  not 

to  create  pitch  or  roll  moments.   The  water  must  act  at  a 

point  x  =  0,  y  =  y   in  the  port  hull,  and  x  =  0,  y  =  -y 

c  c 

in  the  starboard  hull.   Furthermore,  the  amount  of  water 

pumped  is  the  same  for  both  hulls.   The  power  required 

is  the  product  of  equations  (31)  and  (37).   Details  of  the 

calculation  are  presented  in  appendix  E. 

In  roll,  the  control  moments  are  generated  by  adding 

water  to  one  hull,  and  removing  the  same  amount  from  the 

other.   In  order  that  no  pitching  moments  are  generated, 

the  water  must  act  at  x  =  0,  y  =  y  in  the  port  hull,  and 

x  =  0,  y  =  -y   in  the  starboard  hull.   The  power  required 
c 

is  the  product  of  equations  (32)  and  (37 )»  and  is  detailed 
in  appendix  E. 

In  pitch,  the  control  moments  are  generated  by  adding 
water  forward  and  removing  it  aft,  or  vice-versa.   Although 
the  water  can  act  at  any  symmetrical  points  along  the 
x-axis,  it  is  assumed  that  these  points  are  the  hull- 
column  intersections.   An  equal  amount  of  water  thus  acts 
at  each  inter  ection.   The  power  required  is  the  product 
of  equations  (33)  and  (37)  and  is  detailed  in  appendix  E. 

After  the  calculations  in  appendix  E  are  carried  out, 
the  power  required  to  control  the  SPSV  in  heave,  roll,  or 
pitch  is   of  the  form: 

P  =  P  £(BD  cos(wt-T|)  +  #BC  cos(2wt-iy-  X  )  +  #3C  cos(A  -Ty  )] 

(38) 
The  coefficients  are  defined  in  appendix  S.   This  expression 
shows  that  the  total  power  requirements  are  composed  of 
a  constant  component,  a  sinusoidal  component  of  frequency 
w,  ^nd  a  second  harmonic  component  of  frequency  2w. 

The  poxver  required  to  pump  water  out  of  the  tank  is 
positive.   The  negative  power  is  delivered  when  the  water 
re-enters  the  tanks,  but  it  is  assumed  that  this  power  is 
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not  feasibly  recoverable.   As  a  result,  the  average  power 
requirement  over  one  period,  determined  from  equation  (>;), 
is  composed  of  only  the  positive  contribution  of  the  in- 
stantaneous power,  and  is: 

P  SJf  +  §-.  +  K   3C  cos(A-x\)]  (39) 

The  average  pov/er  is  used  to  evaluate  the  feasibility 
of  active  displacement  control  for  a  particular  application. 
Although  peak  pov/er  gives  an  indication  of  the  maximum 
pumping  facilities  required,  and  therefore  capital  invest- 
ment, the  average  power  gives  an  indication  of  the  operat- 
ing costs  over  the  system's  life.   Although  both  peak  and 
average  power  requirements  are  important  considerations, 
only  the  average  power  requirements  will  be  used  for  comp- 
arison in  the  next  section. 
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VII.  DISCUSSION  AMD  CONCLUSIONS 

To  determine  the  feasibility  of  active  displacement  control 
for  a  typical  submarine  pontoon  supported  vessel,  the  power 
required  to  provide  various  amounts  of  control  under  specified 
conditions  is  calculated.   The  characteristics  of  the  SPSV 
are  adapted  from  reference  7  and  are  listed  in  table  1. 

The  first  case  considered  is  that  arising  from  the  desire 
to  have  essentially  zero  motions  in  a  moderate  sea  state.   A 
sea  state  3»  with  a  wave  amplitude  of  7  feet  and  a  period  of 
8  seconds  is  representative  of  this  condition.   The  uncontrolled 
heave  amplitude  is  calculated  from  the  solution  to  equation 
(20-a)  and  is  1  foot,  with  a  zero  phase  angle.   Since  the  motion 
is  to  be  zero,  the  amplitude  factor  (  A  )  is  zero.   The 
resultant  average  power  requirement  is  1930  HP  from  equation 

(39). 

The  next  case  arises  from  the  desire  to  reduce  the  heave 
motion  in  heavy  seas  to  a  level  where  work  can  continue.   A  sea 
state  6,  with  a  wave  amplitude  of  18  feet  and  a  period  of  1^.5 
seconds  is  representative  of  this  condition.   The  uncontrolled 
heave  amplitude  is  9.6  feet.   For  purposes  of  this  discussion, 
it  is  assumed  that  control  costs  and  motion  costs  have  been 
determined  which  result  in  an  amplitude  factor  of  0.35  from 
equation  (27).   In  reality,  the  control  costs  and  motion  costs 
would  be  explicitly  determined,  however  this  determination  is 
beyond  the  scope  of  this  thesis.   With  an  amplitude  factor  of 
0.35,  the  resulting  controlled  motion  is  3.36  feet.   The  average 
power  required  to  maintain  this  heave  amplitude  is  calculated 
from  equation  (39)  and  is  75^0  HP. 

In  the  case  of  a  severe  storn  where  the  wavs  period 
approaches  the  natural  period  of  the  SPSV,  a  resonant  condition 
could  exist.   It  might  be  desirable  to  maintain  the  heave 
motion  within  a  ran<;e  where  limited  work  could  be  performed. 
To  provide  a  comparison  to  the  previous  example,  a  storm 
condition  with  wave  amplitude  of  18  feet  and  a  period  of 
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18.6  seconds,  the  natural  period  of  heave,  was  considered. 
The  resulting  uncontrolled  motion  is  24 .2  feet.   Again  for 
comparison  with  the  previous  example  an  amplitude  factor  is 
taken  as  0.35,  which  gives  a  controlled  heave  amplitude  of 
8.5  feet.   The  average  power  required  to  provide  this  control 
is  28,000  HP. 

If  it  is  not  feasible  to  provide  small  motions  in  the 
resonant  condition,  it  might  still  be  desirable  to  limit  the 
heave  amplitude  to  some  maximum  value.   One  obvious  value  is 
the  v/ave  amplitude.   If  the  SPSV's  heave  amplitude,  in  the  case 
previously  considered,  were  limited  to  the  wave  amplitude, 
18  feet,  the  corresponding  amplitude  factor  is  0.744.   The 
average  power  required  to  maintain  this  maximum  heave  is 
calculated  from  equation  (39)  and  is  9700  HP. 

Table  2  summarizes  the  preceding  examples  and  also  shows 
the  complete  power  equations,  equation  (38). 

The  preceding  calculations  were  based  on  the  ballast 
tanks  being  vented  to  the  atmosphere.   The  power  equations  in 
table  2  show  that  a  significant  portion  of  the  average  power 
arises  from  the  first  term,  or  the  component  with  frequency  w. 
If  the  amplitude  of  this  component  could  be  significantly 
reduced,  the  average  power  requirement  could  be  similarly 
reduced.   Appendix  E  shows  that  the  coefficient  BD  depends 
among  other  things,  upon  the  air  pressure  inside  the  tank.   It 
was  originally  derived  from  a  consideration  of  the  difference 
in  head  across  the  ballast  pump.   If  the  tank  were  maintained 
at  a  constant  pressure  of  14.4  psi,  then  D  =  0  and  the  power 
requirements  are  significantly  reduced. 

The  natural  period  of  the  JPSV  depends  strongly  on  the 
column  diameters.   For  example,  if  the  diameters  were  reduced 
from  25  to  20  feet,  the  natural  period  of  heave  is  increased 
from  18.6  to  23  seconds.   Reduccion  of  the  waterplane  area  by 
reduction  of  the  column  diameters  is  thus  a  means  of  increasi  ng 
the  natural  period  to  avoid  a  resonant  heave  condition.   However, 
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changes  in  the  column  diameters  are  limited  by  other 
considerations.   The  primary  limitation  is  the  "spring 
constant",  or  tons  per  inch  immersion  (TPI)  of  the  3PSV. 
Too  great  a  reduction  in  the  diameters  may  well  result  in 
unacceptable  vertical  displacements  caused  by  addition  or 
removal  of  weight.   This  also  holds  true  for  list  and  trim. 
One  possible  compromise  which  would  permit  small  column 
diameters  at  normal  operating  conditions,  yet  would  provide 
increased  restoring  forces  at  large  drafts,  is  a  column  design 
where  the  diameter  increases  with  submergence. 

To  maintain  buoyancy  with  smaller  columns,  the  hulls  must 
be  enlarged  to  provide  the  necessary  volume.   Structural 
considerations  may  also  limit  the  minimum  column  diameters. 

If  it  is  not  possible  to  increase  the  natural  period, 
active  displacement  control  might  be  considered.   From  table  2 
it  can  be  seen  that  the  average  power  required  to  decrease 
heave  in  a  resonant  condition  is  extremely  large  if  the 
desired  amplitude  reduction  is  large.   However,  if  it  is 
satisfactory  only  that  the  SPSV  not  exceed  the  wave  amplitude, 
the  power  requirements  are  significantly  decreased.   Even  in 
this  situation,  the  power  requirement  is  high,  therefore  it  is 
probabfy  more  economical  to  design  the  SPSV  with  a  high  enough 
natural  period  so  that  resonant  operation  is  avoided. 

The  use  of  active  displacement  control  to  extend  the  use 
of  the  SFSV  into  higher  sea  conditions  appears  practicable, 
although  the  power  requirements  are  still  fairly  high.   The 
greatest  potential  for  active  displacement  control  of  the 
SPSV  is  to  maintain  very  small  motions  in  moderate  sea  condi- 
tions.  The  cost  of  such  a  system  in  terms  of  average  power 
requirements  is  small,  and  although  there  are  practical 
difficulties  in  the  actual  design,  the  concept  should  be 
pursued. 
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TABLE  I 

Characteristics  of  the  SPSV 

Platform  length  185  ft 

Platform  breadth  185  ft 

Hull  length  2k0   ft 

Hull  diameter  25  ft 

Column  diameter  25  ft 

Draft  45  ft 

Clearance  between  water  and  platform               *f*f  ft 

Propulsive  power  6000  HP 
Natural  periods  of  motion 

Heave  18.6  sec 

Fitch  ^2.5  sec 

Roll  55.0  sec 
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TABLE 

2 

Summary 

of  Power 

Requirements 

I 

II 

III 

IV 

A 

7 

18 

18 

18 

T 

8 

14.5 

18.6 

18.6 

lo 

1 

9.6 

24.2 

24.2 

A 

0 

0.35 

0.35 

0.744 

*c 

0 

3.36 

8.5 

18 

p  (pT= 

0) 

1930 

7540 

28000 

9700 

p  (pT= 

14.4) 

68.3 

1150 

15400 

4880 

Condition  I  is  the  state  3  sea. 

Condition  II  is  the  state  6  sea. 

Conditions  III  and  IV  are  the  resonant  conditions. 

The  full  power  equations  are: 

P  =  180  ^32.5  cos(wt  -  270)  +  1.19  cos(2wt  -  270 )J 

P   =  9120  \_   2.2  cos(wt  -  80.5)  +  .123  cos(2wt  -  123.5) 
+  0.0865^ 

P    =  45200  \_0.8?   cos(wt  -  140)  +  0.315  cos(2wt  -  240) 
+  0.24lj 

P   =  45200  \^0.344  cos(wt  -  140)  +  0.103  cos(2wt 


'.0737J 
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DETAILS  Ov    F0RC3  DERIVATION 


The  wave  force  on  a  column  may  be  obtained  from 
ion  (5). 
y-direction  is 


equation  (5).   On  an  element  dz.  ,  the  wave  force  in  the 


Cm  +  A^i  )  \  TTd2 

c      c  op  c  . 
__-t  _ —  ^z. 

m  £y  k            i 


The    total  wave    force   in   the   y-direction   is 
-b 


"  I    \  (mc    + 


Am    )   w2Aek(z:    +  b)sino(sin/r.    dz.  (A-l) 

c  11 


-IH*W) 


The  force  due  to  the  hydrodynaraic  inertia  is: 
-b 

Amc   ( 

— r: —   \  (relative  acceleration  in  the  y-direction)  dz. 

-i.M  u-2) 

The  force  due  to  damping  is: 
-b 

c       r 

-r~   \  (relative  velocity  in  the  y-direction)  dz. 

~<Lb)  x  ^A-5) 

Substituting  the  expressions  for  relative  acceleration 
and  velocity  into  equations  (A-2)  and  (A-3),  integrating 
and  adding,  the  total  force  on  a  column  in  the  y-direction 
becomes: 

Y        =   -\Am    .  >)     .    +  C    .  *1     .1    -    (m        +  2  Am    .  )w  AQ  sin<xeinT. 
ci  \_"    ci    'ci  yi   Lcij  ci  ci  o  i 

-   C    .wAQ  sin*  cost.  (A-*t) 

yi      *o  i 

where       yi.    =    V)+c.i\)-z<D 

xhe  force  on  a  column  in  the  x-direction  is  derived  in 
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a  similar  manner,    and   the   result   is: 

X        =  -  ^Am    .  \         +   0      k     .1    +    (m        +  2  Am    .  )w  AQ  cos^  sinT  . 
ci  L      ci 'ci  xx  5  cxj  ci  ci  o  i 

+  C    .  wA'vi  cos  oncost.  (A-5) 

xx        o  x 

where      fc      .=     t+z0-y.,4-, 
^>  CX  3  c  Jci 

The   derivation   of   the    force    on   a   hull   in   the   y-direc- 
tion  is   similar    to    that    for   a   column.      The   resulting   force    is: 

Y,  .    =  -|Am.  .V},  .    +  H    .n.  .      -    (m,  .    +  2Am..)w  AUsin^sint. 
hi  [       hi  Ui  yx   ^hx]  hx  hx  r  v  x 

-   H    .  wALLsinc\cos  T  .  (A-6) 

yi      r~  a. 

where      y^.    =    ^  +  x^  +  +   z^  cf> 

In  the  x-direction,  the  drag  and  hydrodynamic  inertia 

of  the  hull  are  negligable  compared  to  that  in  the  y-direction. 

Consequently,  the  wave  force  is  the  only  significant  force 

and   is: 

2 

X,  .    =  m,.w   AuLcosoC  sinT  .  (A-7) 

ni  hx         r  i 

The  derivation  of  the  force  on  a  column  in  the  z- 
direction  is  similar  to  that  on  a  hull  in  the  x-direction 
with  the  addition  of  a  restoring  force,  and  is: 

Z        =   -/i-PTTgd2.  t  •    -    (m  w2Q     -   ^TTfr-d2.  )Acost  .  (A-8) 

ex  '     5   cih  c      ^o  \      ex  x 

-he  derivation  of  the  force  on  a  hull  in  the  z-direction 
is  similar  to  that  on  a  column  in  the  x-  and  y-directions, 
and  is: 

-  35  - 


Zhi   ■  -^nhi^'i   +  H,liJ]    -    <<\   +  2Amh)w2ApcosT. 

+  H    .wAUsint.  (A-9) 

zi     n         i 

After  some  simplification,  equations  (A-4)  through 
(A-9)  are  reproduced  as  equations  (17). 
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APPENDIX  B 


J./TAIL3  OF  ^UATIGNS  UT  MOTION 


When  the  force  compoments,  equations  (17),  are  substi- 
tuted into  equations  (18),  the  equations  of  motions  result 
in  the  following  form:  (see  also  reference  1) 


where : 

i 

i 

H  ■ 

3 

9 

* 

and 


The  non-zero  coefficients  are 


~11 
a15 
a22 
a24 

a^4  = 
a,.  —  = 


=  i-i  + 


4m 


4m  z.. 
c  1 

H  +  4m  +  4  Am, 
c       h 

-4m  z.,  +  4  Am.  z, 

c  1       h  h 


=  M  +  4A 


mh 

I   +  4m  z.  +  4Am,z,' 

x     c  2  h  n 

2  ,  a    2 


55 
l66 


Iy  +  4mcz2  +  4Amhx2 

=  I   +  4m  (x2  +  y2)  +  4A 
z     c  c    c 


2 
h~2 


m,  x 


11 


'15 


22 


=  4C 


4z.,C 
1  x 

C  +  H 
x    y 


)_■  =  -4z_C   +  4z,H 
24      1  y     h  y 
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b33  =  S 

V  =  ^CA  +  Vh} 

bKt.  =  *t(C   z?,  +  H  x?J 

55  x  2          z  2 

b~  =  My2C  +  x2C      +  x2H    ) 

ob  ''ex          cy          2   y 


c53  =  ^f  gdc 

ejrt  =  Itf  gd^y2    +  ifmcgz1   -  4mhgzh 

c55   =  ^f  gdcxc    +  H52!   '  Hgzh 


Fn    =   WAS   coso<W-3;^  raw-   Li  m,  w)    sin  wt    +    3   C      cos   wtf 
1  locih  ox  J 

?2   =  WAS   sin*   i  wJ3^  ra      -    UL(m      +  &m    )j     sin  wt   - 

(QC      +  UH    )    cos   wt)  - 

o  y      P  y  J  d2 

F3   =  -WAS  ^Hz  (A  sin  wt    +  [mcwQo   -  *£{  g^-  + 

(mh+2  Amh)w  ulJ  cos   wt? 

Cr                                  2        TT         2  ytandcy   sin*) 

F.     =   WAS   sinoW     -3m  wQn    +    (m  w  Q  -iffgd    )   ~ r= 

4            (.Lcl     c   ^o  *f  »   c  w  sin  c* 

+  (m,  ^Am,  )wU.z    sin  wt  +  (Q,  C   -  z  U^H  )  cos  wtj 

■p  i,    1(,[r   7        -  _       /        2^         TT-,,2.         tan(ky   sin  c<  ) 

F^   a  WAS  Jt  -3mcwi1coso<+    (^w  Qq   -  ^  ^  gdc  ;Xfi  Jc 

w 

+  YphCos«(  +  ^(   [mn2Arnhl    Ks   +   HzKc[!    Sin  Wt 
+      ^CxC06O<  +^e(  ^m^AmJ    Kc    -   RM^)]    cos   wt  J 

F^    =   -WAS  1     3m  w'i  cos*    y      +   Q  C   sinc<   x   tan(kx   coso(  ) 
6  L  L      c      o  Jcoy  c  c 

-  Q  C   cos<*   yctan(kycsino<  )    -  jJ^C  [-^+2  Ara^sin^  KQ 

-  II   sine<    K    )1    sin  wt    +     ,  J>m  wO   sinc<  x   tan(kx   coScO 

ysj  L_c^o  c  c 

+     U.   (    Jn  +2AK,]    sin**    K      +   II   sinc<    K    )J    cos  wtf 
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Written  outt  the  equations  become: 

*n'i +  ai5^'+  hni +  bi5f  =  Fi  (2-2) 

a22T|    +  a2Zf  $+  b^    +  b^d?   =   F2  (B-3) 

*33  ?  +  b33   S  +  °33  S  =  F3  CB-'O 

a2/+  T\   +  %k   ?  +  b24  ^   +  V  ^  +  cMf  ^  =  Fit  (3"5) 

a15  |    +   a^   9  +   b15   £  +   b55  0    +   c55  £  =   F5  (B-6) 

a66*  +  b66*  =  F6  (3"7) 

Since  heave,  pitch,  and  roll  are  of  primary  signif- 
icance, surge  can  be  eliminated  by  combining  equations  (B-2) 
and  (B-6),  and  sway  can  be  eliminated  by  combining  equations 
(B-3)  and  (B-5).   Usin:;  the  notation  D=  —  ,  equations 
(B-2)  and  (3-6)  become: 

D(a--D  +  b_.)fc   +  D(a-_D  +  bnc.)0  =  F-  (3-3) 

11     11  *      15     15      1 

D(a15D  +  b15)  I    +  (a55D2  +  b^D  +  e^)  6   =  F^  (B-9) 

Solving  simultaneously  for  &  : 

[(allD  +  bll)(a55D2  +  b55D  +  c55]   "  D(a15D  +  b15)2]  ^  = 

(allD  +  bll)?5  "  (a15D  +  b15)Fl        (B-10) 

Performing  the  required  operations,  and  noting  that: 

F.  =  f.  sin  wt  +  f.  cos  wt 
i    is  ic 

where  the  coefficients  are  given  above,  equation  (3-10) 

becomes: 

•»•      ••       • 

C13  ^+  G12  ®  +  Cll  ^  +  G10  ^  =  Jlssin  Wt  +  JlcCOS  wt 

(B-ll) 
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where : 

2 
C13  =  alla55  "  a15 

C12  =  allb55  +  blla55  "  2a15b15 

Cll  =  allc55  +  bllb55  "  b15 

C10  =  bllc55 

Jls  =  "allwf5c  +  bllf5s  +  a15wflc  "  b15fls 

J-   =  annwfc   +  b-.-f..  -  a..  cwf .,   -  bnc.f. 
lc     11  5s    11  5c    15  Is    15  lc 

Noting  that  equations  (B-3)  and  (B-5)  are  similar  to 
equations  (B-2)  and  (B-6),  the  differential  equation  for 
roll  can  be  written  by  inspection: 

•••  ••  • 

G23  $+  C22  CP+   C21  <P  +   C20  *?  =   J2sSin  Wt    +   J2cCOS  wt 

(B-12) 

where : 

r  2 

23   =  a22aMf  "  a24 

C22  =  a22b/fif   +  b22a42f  -  2a24b2^ 

2 
C21   =   a22C^   +  b22bM+   "   b2*t 

C20   =  b22CMf 

J2s   =   -a22wf^c    +  h22hs    +  a2^fWf2c   "   b2<ff2s 

J2c  =   a22wf^s  +  h22fkc   '   a24wf2s  "  b24f2c 

Summarizing,  the  equations  of  motion  for  the  SPSV  in  heave, 
pitch,  and  roll  are: 


a,,  *s  +  b_,  \   +  c-,..  K    =  f,  sin  wt  +  f_  cos  wt    (B-13) 
33      33      3s         3c 


33 

C.,  0  +  C.0  Q  +   Cnn  £  +  C. _  g  =  J,  sin  wt  +  J_  cos  wt    (B-l4) 

13      1^      11  10  Is          lc 

0       <P*  +  C22  cp  +  C21  <$>  +  C20  ^  =  J2sSin  Wt  +  J2cC0S  Wt    ^B-15) 
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APPENDIX   C 
DETAILS   OF    SOLUTION   OF   E^UATIOIIS   OF   I-'.OTILN 

The  equations  of  motion  in* heave,  pitch,  and  roll, 
equations  (3-13),  (3-1^),  and  (3-15),  can  be  solved  by  assum- 
ing solutions  of  the  form: 

$    =  I      cos(wt   -  ls)  (J-l) 


9=9    cos(wt   -  t9) 


(C-2) 


Cp  =   cp     cos(wt    -  ?4>) 


(C-3) 


Inserting  these  relationships  into  the  equations  of  motion, 
equations  (20),  results  in  expressions  which  determine  the 
magnitudes  and  phase  angles  of  the  3PSV  motions. 
Substituting  equation  (C-l)  into  (20-a): 

3  U-a  w  +  c   )  cos(wt  -  Zj)  +  (-b  w)sin(wt-  £^)J  = 


f,  sin  wt  +  £-,   cos  wt 
ys  3c 


(C-4) 


Combining  sine  and  cosine  terms  on  either  side  of  equation 
(C-k )    results  in: 


"$  R^  cos(wt-  |5^  )  =  F  *cos(wt-  ^) 


(C-5) 


where : 


K-  =  K-a-^w^  +  c„ ) 


2       N2 
;33 


"b33w)2J 


>. 


h=  ^+  t 


an 


-U    "V* 


-aw   +  c 
33      j? 


I  3c. 
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7or  equation  (C-2)  to  be  true  for  arbitrary  values  of  time: 


and 


■S 


h- 


W 


R 


(J-6a) 


ft- 


or 


r 


tan 


-1 


-b^w 


] 


-a^w   +  c33 


(C-6b) 


aquations  (C-2)  and  (C-3)  are  treated  in  a  similar 
fashion,  the  presence  of  a  third  derivative  adding  one 
additional  term.   Substituting  equations  (C-2)  and  (C-3) 
into  (B-l*+)  and  (3-15),  simplifying,  combining  sine  and 
cosine  terms,  and  solving  for  the  amplitudes  and  phase 
angles,  the  resulting  expressions  for  pitch  and  roll  become 


o    ^ 


*e=    ft-  * 


an 


-1 


3      1 

C13W  -  °11W 
-C12w2  +  C 


where : 

R 


6 


-   [(-C12"; 


10 


+  (C,,w3  ■  C,,w)2j 


)4 


13 


11 


2  J4 


F*  =  U-,   +  J,  ) 


0 


Is 


lc 


(C-7a) 
(C-7b) 


ta  =  tan    T- 


9 


lc 


and  for  roll 


To    R^ 


an 


C23W  "  °21W] 

2 
'C22W   +  C20 


(C-tta) 
(C-Sb) 


-  k2   - 


where 


R4  -    [(-C22"2    +   °20)2    +    (C23w3   "   C21W)2] 


a 


F«  =    (J2s*J2=>" 


r 


V«>  = 


=   tan 


-1 


is 
r2c 


.43- 


APPEKJIX  D 
DETEl^'iINATICN  OF  COilTROLLED  H0TI0H5 

To  determine  the  controlled  motions  of  the  SPSV,  the 
control  force  is  added  to  the  exciting  force  on  the  right 
hand  side  of  the  equations  of  motion,  equations  (20)  as 
illustrated  below  for  heave. 


a,.,  X  +  b,_  \  +   c-,v  \  =  f,  sin  wt  +  f,  cos  wt 


33 


33 


33 


3s 


3c 


U  ^  cos(wt  -  S* ) 


(D-l) 


where:   U..  -   (  S       -    S     )  from  equation  (26), 

>  r>  J  o    J  c 


The  controlled  motion,   S  =  S   cos(wt  -  i.  <  )  is 

c  >c 

substituted  into  equation  (D-l)  resulting  in: 


1    R^cos(wt   -   (i     )   =  T  ^    cos(wt   -  K-) 


(D-2) 


where : 

R 


S    =  [(-a35w2   +  c53)2   +   (-b33w)2J 


)i 


01.   -    £S=   +  * 


an 


-II 


-b33W 


-a^w      +   c33  A 


T  ^    =  Tf^  +   U^  +   2FSU^   cos(    )k  -    Ss)j 


/s 


-1 


K,   =  tan        1F     c 


F  ^sin  ys     +  u- 


sx 


n    fc. 


os    ^    +  U^  sin    ^ 


and  F ^  and  £*  are  defined  in  appendix  C.   Equation 
(D-2)  results  in: 


s--*t 


and  Z .      -    K  <.  -  tan 

>C       a 


-1 


"b33W 


-a33w  +  c33 


(D-3a) 
(D-3b) 
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Expanding  equation  (D-^a)  and  squaring: 

2   «?   f^  +  u2  ,  zf,u  .«(>,  -S  ) 

A  A 

From  equation  (C-6a),  "$   =  — —  ,  therefore 

O    "5 

2      2         U-5  A 

Substituting  U*  from  equation  (26)  into  equation  (D-5), 


2    „  3  o  q3   ,  r2   r2^ 


Define  an  amplitude  factor  (A)  as  the   ratio  of  the  control- 
led motion  to  the  uncontrolled  motion: 

/\  5  -£■£  (D-7) 


I 


o 


'/hen  equation  (D-7)  is  substituted  into  equation  (D-6)  and 
the  expression  simplified,  the  result  is: 


2         R"S  PS     ,.    a2^ 


/\^  -  2  >     (1  .  /\£)*  cos(  ^  _  J  )  -  1  .  0 

5  /   n    \ 


1  +l"  (D-6) 


The  motions  of  pitch  and  roll  are  treated  in  a  similar 
manner,  with  the  subscript  "S  replaced  by  6*  or  op  as 
applicable.   The  resulting  equations  are  similar  to  equation 
(D-5).   Without  subscripts,  Lhe  general  equation  can  be 

-  *f5  - 


written: 


sr  -  k  (i  - 


A2r  coS(  ^  -  s> ) 


-1    =   0 


(D-9) 


where: 


K   =   2 


-2lL 


w 


Equation  (D-9)  is  a  general  expression  relating 
controlled  motions  of  the  SP37,  resulting  from  the  optim- 
al control  force,  to  the  phase  angle  of  the  control  force, 
for  given  values  of  the  control  cost  and  motion  cost. 


-  kS   - 
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DETAILS  F  PO'itESi   CALCULATIONS 


-he  po^/er  requirement  in  heave  is  the  product  of  the 

heave  :nass  flow  rate,  equation  (3D,  and  the  head,  equation 

(37).   Since  heave  is  the  only  motion  considered,  @  =(0=0, 

The  mass  flow  r,j.te  for  each  hull  is: 

i-  =  -^PTTwd2(  \     sin  wt-  £<  )  -  $  sin(wt-  ?   ))      (E-l) 
■*     '     C  "  o         »       c  >c 

The  head  seen  by  the  port  hull  is: 

PT  -kz 

H  =  Z      -   -—   +   Ae    ""  °cos(wt    +   ky   sine*  )    -     X    cos(wt-  £.    ) 
p        o        fg  Jc  Jc  Sc 

d2 

-57 -\Z   cos(wt-f-)    -    X    cos(wt-  *       )I  (E-2) 

T  "S     *- 

The  head  seen  by  the  starboard  hull  differs  only  in  the  sign  of 
y  a   -he  power  requirement  is  therefore: 

P  t  =  mH  +  mH  =  m(H  +  K  )  (E-3) 

I  J  p    *  s    s  P    s 

./hen   the    expressions    for    head   and  mass    flow   rate    are    substituted 
into    equation    (E-3)>    and    the    equation   is   expanded    and    simp- 
lified,   the   power   equation   becomes: 

r  3*cs 

P.   =   r       \  3j  IK   cos(wt-YW  )    +  — r —  cos(2.;t-  r\    -   A-^) 

:os(    \s-    ^  )J  (E-*f) 


~5^S 


where  : 

2 

0"i  '  CO 


3  ,    =[l    ♦     A,2   _   g  Ascos(^c    -  £sf}* 
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APPENDIX  E 


DETAILS  OF  POWEH  CALCULATIONS 


i'he  power  requirement  in  heave  is  the  product  of  the 
heave  mass  flow  rate,  equation  (3»Dj  and  the  head,  equation 
(37).   Since  heave  is  the  only  motion  considered,  0  =(0=0. 

C    1  c 

The  mass   flow  rate    for   each  hull   is: 

nus  -/2pirwd2(  }    sin(wt-E-)   -    I    sin(wt-  ?      ))  (E-l) 

*  •  c     *  o  »  c  "Sc 

The    head   seen   by    the   port    hull   is: 

T  — kz 

H   =   z      -   -—  +   Ae        °cos(wt    +   ky   sino<  )    -     X    cos(wt-  £-    ) 
p        o        fs  c  c  >c 

d2 

-^ \l    cos(wt-f-)    -     X    cosiwt-£t     )1  (E-2) 

T"S     *- 

The   head   seen  by   the    starboard   hull   differs   only   in    the    sign  of 

y    .      i'he   power   requirement    is    therefore: 
c 

P  t    =  mH      +  mil     s  m(H      +   H    )  (E-3) 

S  tp  *  s  SP  s 

.v'hen  the  expressions  for  head  and  mass  flow  rate  are  substituted 
into  equation  (E-3)»  S-^d  the  equation  is  expanded  and  simp- 
lified, the  power  equation  becomes: 

r  3ics 

F   =  P   \  Bj  D<  cos(wt-Y^)  +  -5 —  Cos(2wt- Y^-  h^) 

+   -|-i  cos(  Xs-  t^  CE-4) 

where : 

p  =evra2vilz 
o-\      »   c  •*  o 


s  -Jl  +  f\l   -  2  Ascos(^c  -C*fl* 


-  <*7  - 


n*- 


=  tan 


-1 


A«cos    f          -   cos    Z 
i Lc ± 


\_Assin  &sc  "  sin  ^\ 


v-fc<-.-ri> 


A      -kz. 


TTd. 


cos(ky   sin  ok  )    +    (5-7 —  -   1)  /\.cob  I 


2  A 


TS 


2ATS  SJ  ^   2AT^  \\  Sc 


■Sc 


TTd 


2., 


sin  £ 


TS 


A<  =   -    tan 


TTd 

"2h 


2 

£ 


(^T1  -  1)  Assin£ 


TTdc 

sc  -  2a77  sin  Z  S 


A     -kz 
Jo 


TTd 


TTd; 


-  e~    Z"cos(ky   sinc<  )    +    (  -   1)  A*cos  £  f      - 


air:  cos  ei 


The  power   required   to   control  roll   is   determined   from 
the   product    of   equations    (32)    and    (37).      The   mass    flow  rate 
for   each   hull    is: 

m.=  -#fTTwd  y   V<p  sin(wt-I^)   -  d^  sin(wt-  £.    )l  (£-5) 

The  head  seen  by  the  port  hull  is: 

Pip  — kz 

H      =   z     -  •=—  +  Ae        °cos(wt    +  ky   sin  <* )    +  y  <P  cos(wt-  ?«   ; 
p  o        fg  Jc  ^cTc  Vc 

TT  rl 

"   2lT/c  [l>0=o=(",t-£?)    -  ^cosfwt-^j  (E-6) 

The  head  seen  by  the  starboard  hull  differs  only  in  the  sign 

of  y  in  the  third  term.   The  power  requirement  is  therefore: 
c 


P4>=m<p(HP  +  V 


(E-7) 


-  if8  - 


Substituting   the    expressions    Tor   head   and   raass    flow   rate, 
t^e   resulting  power   equation   becomes: 

+  >23^C^ccs(  >^-  \^)]  (E-o) 


where : 


P   A=PTfd  w(v    <£    ) 
o$    i        c      *  c  ^o 

A^cos  2^>c    -    C03  2^"" 


A« 


=    tan 


A^sin  £^c    -   sin  €^ 


1        r  ^ 

yc  ef0      o     f  g 


c     =  J 


-kz 


ird. 


:   To 


e     "  °cos(ky    sin  <x  )    +    (— +   1)  A<dc°s  2- 


<9 


TTd' 

2i^:cos  2* 


«? 


*n  2  T    Tt  d^ 

♦J       +     [(^+l)Vin 


TT  d2 
c 


% 


'<?c        2Am 


sin  £, 


<? 


A<d=  -  tan 


-1 


2  2 

TT  d  TT  d 

( ~ ,         +   D/V^sinf         -  sin   £4) 


.  Trd^ 

— V-  e        °cos(ky    sin  *  )    +    (^-r-2-   +   l)A$cos£ 
yc«9o  C  2>1T4>  Y 


TTd 


4>c    2A 


co.;  £ 


T  <? 


4> 


The  power  requirements  for  pitch  are  determined  in  a 
similar  manner,  and  the  power  equation  is: 

Pg  =  P   pgD^ccsCwt-T^g  )  +  ;l3^J^cos(2wt--H  -  Ag) 

+  ;iBdctfcos(/e-y]  (e-9) 


-  ^9  - 


where : 


P     _  =fTTd2w(x  Q  )2 

Q&        '  C  CO 


^9 


=  [l   +  hi   -   2A6>cos(E6>c    -^ 


ajfi 


^9  =  tan 


-1 


Agcos  £         -   cos  £, 


0c  ""  ^O 


A*sin  lec  -  Gin  L°  . 


9 


e 


i 

c  ^o 


Pm 

(z      -    o1) 


\ 


•kz, 


TT  d 


x     0 

c  ^o 


(SI 


-*  +  dA 


TTdc  *-%2 


f    'c  A  ^c  "I 


T0 

/2 


A„  = 


0 


-tan 


-1 


^d2 
(T +   1)  Apsin  I 


TTd 


hk 


T  9 


ec  "  5a 


>in  £ 


T0 


0 


M.o  -r.z0 


l-  c  w  o 


2  2 

TTd  TTd 

+    (j-r-S  +  1)  A^cos  t        -  -^ 


cos  2, 
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